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We characterize the motion of charged as well as neutral tracers, in an electrolyte embedded in a 
varying section channel. We exploit a set of systematic approximations that allows us to simplify the 
problem, yet capturing the essential of the interplay between the geometrical confinement provided 
by the corrugated channel walls and the electrolyte properties. Our simplified approach allows us to 
characterize the transport properties of corrugated channels when a net flux of tracers is obtained 
by keeping the extrema of the channel at different chemical potentials. For highly diluted tracer 
suspensions, we have characterized tracers currents and we have estimated the net electric current 
which occurs when both positively and negatively charged tracers are considered. 


I Introduction 

Both in biological situations and in synthetic devices charged particles are transported across a channel whose 
radial section ranges from the micro- to the nano-metric size. In such systems the confining walls are usually charged 
and the electrolyte organizes around the walls to screen their charges. As a result, an inhomogeneous electrostatic 
field develops inside the channels that affects the transport of charged tracers along the channel, leading to a charge- 
dependent channel permeability. Until now, the control of such currents has generally relied on applying an external 
driving whose origin is usually either of electrical or hydrostatic origin. Tuning the external drive allows, for example, 
for the control of particle currents as it happens in sodium-potassium pumping in neurons [l[. However, the shape of 
the geometrical confinement provides, in itself, an alternative route to transport and current control. Different groups 
have shown that, in the presence of external forces, the local variation in channel section can induce novel dynamical 
regimes such as particle separation , cooperative rectification and negative mobility just to mention a few 
among others [6l-l8|. In particular, it has been shown that the inhomogeneous distribution of particles along the radial 
direction Q can modulate the overall current as it happens for neutral tracers under an external field, such as gravity. 
Moreover, when an electrolyte is embedded in a varying-section channel, novel dynamical regimes, such as particle 
separation and negative mobility, can be attained when the Debye length (characterizing the decay of the electrostatic 
potential from the channel walls) and the channel section are comparable 0. Such a regime differs from previously 
studied scenarios SEl when the Debye length was much thinner than the channel’s cross section. 

Membrane permeability depends on the membrane structure and also on how molecules diffuse through it. It is 
in general not obvious how to disentangle to what extent the measured permeability is essentially controlled by the 
morphology or by the dynamics. Therefore, we will not consider the permeability associated to the flow through a 
pore in response to applied pressure gradients, as is usually analyzed in the context of porous materials. Rather, in 
this piece of work, we will study the permeability of a corrugated channel driven by a difference in chemical potential 
at the ends of the channel. We will focus on the interplay between diffusion and local changes in the pore geometry, 
and will concentrate on the effect that channel shape has on the diffusion of ions inside the channel. In particular, we 
will not consider the permeability associated to the flow through a pore in response to applied pressure gradients, as 
is usually analyzed in the context of porous materials. Rather, in this piece of work, we will study the permeability 
of a corrugated channel driven by a difference in chemical potential at the ends of the channel. Such a regime has the 
advantage that the fluid ions diffuse through can be considered as at rest therefore simplifying the overall dynamics. 

The permeability in ion channels, and in membranes, quantifies permeation across permeable substrates. In par¬ 
ticular, the permeability coefficient is an involved function that depends both on intrinsic properties of the soluble 
molecules, such as their concentration and diffusion coefficient, on properties of the channel, such as its size and geom¬ 
etry, and in the interaction between the molecules and the channels, quantified for example through the partitioning 
coefficient (which quantifies the affinity of the molecule to the channel) or a slip coefficient. The Goldman-Hodgkin- 
Katz theory accounts for such ingredients and the electrostatic interactions between molecules and charged ions and 
predicts a strong, non-linear dependence of the permeability when a voltage difference is applied through the channel. 
This theory emphasizes the relevant contribution of the motion of the ions inside the channel to the channel perme¬ 
ability The results we will identify for the effective tracer diffusivity along the channel will provide fundamental 
understanding on the impact that channel shape has on charged ion permeabilities in both channels and membranes. 
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FIG. 1: Electrostatic field inside a varying-section channel whose bottleneck half-amplitude ho — h\ is comparable with the 
Debye double layer thickness, A. 


Specifically, we shall show that the shape of the geometrical confinement provided by the channel walls as well as the 
electrolyte properties, captured by the Debye length, lead to a significant modulation of the overall channel perme¬ 
ability for both charged and neutral tracers. In particular, we find that the relative position of the bottlenecks as 
compared to the reservoirs at the end of the channel plays a key role in the channel permeability, therefore providing 
additional ways to tune the net flux along the channel. The structure of the text is as follows: In Section II we will 
derive the Fick-Jacobs equation for charged tracers moving in a varying-section channel, in Section III we will present 
our results and finally in Section IV we will provide our conclusions. 


II. Electrolyte embedded in a corrugated channel: a Fick-Jacobs approach 


We will consider a symmetric, z — z electrolyte embedded in a negatively charged channel of varying half-section 
amplitude, h(x), although analogous results can be obtained for an insulating channel, subject to a gradient of tracers’ 
concentration. In order to describe the dynamics of suspended charged tracers, one needs to analyze the diffusion 
motion of the tracers and the electrostatic potential inside the channel, which constitutes a formidable task. However, 
the analysis can be significantly simplified for highly diluted ionic concentrations and small C potential on the channel 
walls, i.e. /3^e <IC 1 where = kgT is the inverse temperature, with kg the Boltzmann constant and e the 
elementary charge. In such a regime, the Poisson-Boltzmann equation, which determines the electrostatic potential in 
thermodynamic equilibrium, can be linearized and simplified to the Debye-Huckel equation for which the electrostatic 
potential decays exponentially with the distance from the wall over the characteristic length scale, A = k~^, namely 
the Debye length. We will also consider the regime where the charged solute suspension is highly dilute (tracer limit), 
when its dynamics is described by an advection-diffusion equation. In order to gain insight in the properties of the 
channel permeability upon variation of the channel geometry, we will assume that the channel half-section, h{x), 
varies slowly, i.e. dxh{x) <g; 1. Such an assumption allows for a projection of the convection-diffusion equation, which 
determines the dynamics of tracers inside the channel, to an effective \D equation, where the varying-section of the 
channel enters as an entropic effective potential. This approximation, called Fick-Jacobs, has been used [ 1 , 13, 14| 
and validated [^, 15-17| in many different scenarios. 


We will consider a channel, of length L, whose section varies only along the a;-direction and it is constant along the 
z coordinate, with a simple, periodic shape whose half section is determined by 


h{x) = hoi + hi cos 


/ 2'kx 


( 1 ) 


where /iq and hi correspond to the average and maximum channel modulation, respectively. In turn, (j) determines 
the relative position of the channel bottleneck with respect to the connection of the channel to the reservoirs it is in 
contact with. For <() = 0, at x = L/2 the channel shows a bottleneck while for (p = tt the minimal aperture of the 
channel occurs at x = 0. 
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The dynamics of tracers embedded in an electrolyte, quantified by its probability distribution inside the channel 
Pa{x,y, z,t), with a = ±,0 denoting positive, negative and neutral tracers, is governed by the convection-diffusion 
equation, which in the overdamped regime, reads 

dtPaix,y,z,t) = npV ■ [Pa{x,y,z,t)VUa{x,y,z)] + DV^Pa{x,y, z,t) (2) 


where D is the tracer diffusion coefficient and Ua is the total potential acting on the tracers, 

, „ ^'1 _ / Qaeipix, y, z), \y\ < h{x) & |z| < L^/2 

“ \ oo, \y\> h{x)oi\z\> Lz/2 

that is periodic along the longitudinal direction, Ua{x,y,z) = Ua{x -\- L,y,z), and confines the particles inside the 
channel, ijj stands for the electrostatic potential inside the channel that, due to the high dilution of tracers is supposed 
to be unaffected by them, and qa corresponds to the tracer valency while e is the electron charge. We assume that 
all tracers have the same intrinsic diffusivity, D. Hence, the differences in tracer diffusion we will describe later will 
then be essentially due to the environment in which the tracers will move. The electrostatic potential, obeys the 
Poisson equation 


dli’ix, y, z) + dyipix, y, z) + dlp^{x, y, z) = ( 3 ) 

where e is the liquid dielectric constant and the value of the electrostatic potential at the interfaces depend on the 
conducting nature of the channel walls. In the Poisson equation pq{x,y,z) stands for the local charge density of the 
electrolyte. At equilibrium, the charge density is derived using the Boltzmann distribution of ionic densities inside 
the channel. When the channel walls are smoothly-varying, dxh{x) 1, and assuming lubrication dx'tp(x,y, z) <C 
dyip{x,y, z) and d'^'ip{x,y, z) <C dyip{x,y, z), we can reduce the Poisson equation to an ordinary differential equation 
for the potential ^|J along the channel. At a solid, homogeneous [l^ interface, the electric field is perpendicular to the 
solid wall. As a result, for a varying-section channel the effective field along the channel corresponds to the projection 
of the electric field along the wall. However, for a smoothly-varying channel amplitude, the projected electrostatic 
field reads: E = Eocos{0), with 9 = arctan the local channel slope. Since, dxh{x) 1, the corrections on 

the electric field due to changes in the channel section are of second order in dxh{x) and can be safely neglected in 
the following. For low salt concentrations and small C potential on channel walls, we can further simplify the Poisson 
equation by linearizing the charge density pq(x,y) ~ po (1 “ (ize'tp{x,y)), hence getting 


'i/j{x,y,z) 


^ cosh(fcy) 

^ cosh{kh{x)) 


( 4 ) 


for a channel made by conducting walls or 


ilj{x,y,z) 


a cosh(fcj/) 
2e sinh(fcft,(a::)) 


( 5 ) 


for an insulating channel characterized by a constant surface-density of electric charge a. Such an assumption, 
known as Debye-Hiickel approximation, allows to identify the screening length, k~^, of the electrostatic potential as 
= jSzepa/e where po is the number density of the ions in the electrolyte and z their valency. The approximation 
made for the electrostatic field reflects in the Debye length, A, that results constant up to second order in dxh{x). 
Finally, for dxh{x) 1 we can approximate the transverse profile of the probability distribution function (pdf), 
Pc,{x, y, z, t), of a tracer of valency Qa by its profile at equilibrium, i.e., we can factorize the pdf by assuming [20|: 


Pa{x,y,z,t) 

/3Aaix) 


^-0qa,eip[x,y) 

^ l'Lz/2 

— In —- / 

2hoLz J-L^/2 





( 6 ) 

( 7 ) 


After integrating over the channel cross section we arrive at 


(ipa{x,t) 


dAa{x) 


^xPa ix^ ^') 




dx 


( 8 ) 
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This expression encodes both the confining as well as the electrostatic potential in the effective potential A{x) whose 
shape, see Eq. [3 resembles that of an equilibrium free energy. Since all the quantities of interest are independent of 
z, without loss of generality we can assume /_2 /2 dz = 1 and consider all quantities per unit of transverse length, 
Lz. Defining the average, a;-dependent, electrostatic energy as 

(Va^ix)) = (9) 

J — h{x) 

from the definition of Aa{x) we can define the, dimensionless, entropy along the channel as ksTSoiix) = {Va{x)) — 
Aa (x ), from which we can define 


Saix) = In 


1 rHA 

2«o J-h{x) 


- 13 {Va{x)) . 


In the linear regime, I3qaeijj{x,y) <C 1, we can linearize the last expression getting 

'2h{x )' 


S{x) ~ In 


2hn 


( 10 ) 


( 11 ) 


where the entropy has a clear geometric interpretation, being the logarithm of the space, 2h{x), accessible to the 
center of mass of a point-like tracer. Accordingly, we introduce the entropy barrier, AS", defined as 

AA = ln(^), (12) 

which represents the difference in the entropy evaluated at the maximum, h^axi and minimum, hj^im of the channel 
aperture. 


Ill Tracer dynamics 


A net, constant flux of tracers, of magnitude J, characterizes the steady state motion of tracers in the channel 
whenever there is a difference in the chemical potential of the two baths the channel is in contact with. Since there 
is no electrostatic potential difference between the reservoirs at the two channel ends, and the tracer density is small, 
the tracer chemical potential in the reservoirs reads fj,a = ksT In pa- li Pa differs at the ends of the channel, the 
associated A/io, will lead to a tracer flux along the channel. We are interested in understanding the dependence of 
channel permeability on the channel geometry as a result of the imposed chemical potential difference. 

From Eq. ([8]), in steady state we can express the tracer density per unit length profile inside the channel, na(x), as 


,{x) = 


— p-3A^(x) 




(13) 


where we have used that for the tracers Pa{x) and na{x) are proportional to each other [2l|. From na{x), the local 
average tracer density can be finally obtained as pa (x) = 2 h(xj > relation provides the natural link with the tracer 
densities in the reservoirs. Specifically, we will consider that the reservoirs at the two ends of the channel are kept 
at tracer densities pap and Pa, 2 , corresponding to an imposed chemical potential difference Apa = Imposing 

Ha (—j) = 2pa,ih(—j) and ria (-j) = 2paph{^) we determine the two constants in Eq. (fT^ 'l 


n. 


= 2h{L/2)p,,pe^M-^) 


(14) 


Ja 


-2Dh{L/2) 


{Pa,2 — Pa,l) 






(15) 


where in the last expression we have exploited the symmetry of the channel, h{—^) = h{^) that, using Eq.([3), 
implies Aq,(—- j) = Aa{^). Eq. ITSl provides the dependence of the tracers’ mass flow due to both the gradient in 
tracer concentration, Vpa = {Pa ,2 — Pa,i)lL, and the geometry of the channel. We can rewrite Eo. ITKlas: 


Ja = -2Dh {L/2) \7paXa 


( 16 ) 











5 


where the geometrical dependence on the flux is encoded in the dimensionless parameter 




Xa = 


j. ^pAa(x)fl^ 


( 17 ) 


According to Eos. 1161171 Xa captures the dependence of tracers’ flux upon both channel geometry and tracers’ charge. 
Moreover, Eg. ITTI shows that Xa is independent on the magnitude of the drive, Vpa, responsible for the onset of the 
flux. In particular, when Xq > 1, channel corrugation enhances tracers’ flow as compared to the case of a flat channel 
hence leading to a larger channel permeability while a reduced permeability will be obtained for Xq, < 1. 

Eq. (HZD shows that for constant channel sections Xa = 1 irrespectively to tracers’ charge. Therefore, according 
to Eq. [16] when the channel section is constant tracers’ flux, and hence channel permeability, are insensitive to 
tracer charge, provided they undergo the same chemical potential gradient, Vp^- Tracer charge will determine the 
partitioning of the tracer inside the channel and its distribution in the channel section, but will not affect its effective 
diffusivity. 

The dependence of Xa on the channel varying geometry is quite involved, as can be appreciated in Eq. [T71 We can 
gain insight into the impact that channel corrugation has on its permeability by analyzing the particular case where the 
channel amplitude varies linearly, namely h{x) = ho — hi\x\/L. In this geometry, j5A{x) ~ — \n.{ho/L) + {hi/ho)\x\/L. 
This regime can be achieved, for example, if the channel corrugation is weak and the electrostatic potential varies 
linearly. Or, more generically, one can consider a channel whose section does not change strongly and hence one can 
identify the coupled effect of channel corrugation and electrostatic potential variation through an effective, uniform 
variation along the channel. 

For neutral tracers this regime corresponds, exactly, to a channel with a linearly increasing cross section. For 
charged tracers, in the regime of a very narrow double layer, kho 'S> 1, tracers are essentially uniformly distributed in 
the channel section; accordingly, they will behave as neutral tracers. Under these assumptions we have 


ehMAdx = 2 / = 


g/3hi/(2ho) _ \ 


(18) 


Substituting the last expression in Eq.jTTjand using the fact that A (jJ = A (—j) = In/iQ + hi/{2ho) we get 

1 


Xa — 


1 _ g—/3/ii/(2^o) 


(19) 


Finally, recalling that in this case ^AA = hi/{2ho) the last expression reads 

“ 1 _ g-/3AA 

Interestingly the last expression shows that the permeability does not depend symmetrically on AA. Hence, the 
permeability of a concave channel, where the minimal aperture is in the center of the channel, will be smaller than 
the complementary channel where the center of the channel has the largest aperture. For a channel with the same 
mean section, the difference in the permeability for channels with opposed curvature will differ exponentially with the 
channel corrugation. 

It is insightful to analyze the limits 


lim Xa 

= 1 

(21) 

AA^O 

. lim Xa 

= BAA 

(22) 

AA—¥oo 

A 1™ Xa 
AA^ — oo 

= 0, 

(23) 


which clearly show the strong dependence of the permeability on the effective, free energy barrier AA. This free 
energy barrier, which encodes the change in channel section and the local modifications this induces in the electrostatic 
potential, controls the impact of the environment on the net transport properties of charged tracers. 


IV Results 


In order to characterize quantitatively the dependence of the channel diffusivity as a function of tracer charge 
and for different channel geometries and electrolyte concentrations, we have numerically evaluated Eq. 1171 For the 
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FIG. 2: Tracers permeability, a: as a function of the phase (j) for positive (red) and negative (blue) tracers, normalized by the 

current of neutral tracers, for different values of the corrugation, AS = 1.1, 2, 3.7 for squares circles and triangles respectively 
with kho = 10. b: Xa for charged tracers, normalized by the permeability of neutral tracers xo, as a function of the Debye 
length, A, expressed in units of kho for different channel geometries: (f) = 0 {(j) = 0.5) is characterized by circles (triangles) 
while the color code captures the charge of the tracers red (blue) stand for positive (negative) tracers while marker size is 
proportional to the amplitude of the corrugation AS = 2, 3.7. c-d: electric current, as defined in Eq. 1241 calculated from the 
data of the corresponding above panels. 


simple channel profile, characterized by Eq. o, neutral tracer diffusion does not depend on channel shape. Hence, 
the deviations we will analyze emerge from a correlation between channel corrugation and the local variations in 
electrostatic potential. As shown in Fig. [2ja, Xa depends on channel shape and the relative position of the maximum 
aperture with respect to the channel ends, quantified by the phase shift (j). The variation indicates that Xa, and hence 
the permeability, depends symmetrically on the distance of the channel maximum with respect to the channel ends. 
The impact of channel shape on tracers’ dynamics has opposite trends depending on tracers’ charge. Negative tracers 
experience a net repulsion form the solid surfaces and will preferentially accumulate in the channel center. They 
display a maximum in permeability when the bottleneck is at channel edges, while the minimum in the permeability 
happens in the opposite situation. Positive tracers exhibit a complementary behavior. These tracers accumulate 
preferentially at the channel walls and benefit from the shape when the charge accumulated at the two walls do 
not interfere with each other. The non-monotonous behavior of the permeabilities leads to different scenarios where 
the ratio between the permeabilities of charged and neutral tracers can vary significantly. We can identify different 
channel configurations for which one of the three species has the largest permeability. For example, positive tracers 
(hence attracted from the negatively charged walls) experience the maximum permeability when the channel shows 
a bottleneck at a; = 0, while negatively charged tracers are fastest in the case of bottlenecks at channel edges. 
Alternatively, for cj) ~ 0.4,0.6 neutral tracers are faster than both positively and negatively charged tracers. 

Charged tracer permeabilities are strongly affected by the electrolyte properties, which are captured by the Debye 
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length, A, and depend on the combined properties of the channel and the electrolyte, as it has been shown for the 
case of electrostatically driven tracers Q. Fig. [2]b shows that if A and the channel average amplitude, are not 
comparable, the modulation in the permeability due to the geometrical confinement vanishes and both positively 
and negatively charged tracers experience the same permeability. For fc/iQ <C 1 the electrolyte is uniform in the 
channel section, while for kho ^ 1 the electrolyte is localized at the wall and the fluid inside the channel is essentially 
neutral. In both cases all tracers behave analogously and become insensitive to channel corrugation. In contrast, 
when kho ~ I the permeabilities shows maxima (minima) according to tracers’ charge. The regime fc/iQ ~ I can 
be obtained in nanofluidic devices for which the Debye length, that is typically A ~ 1 — 10 nm) matches the typical 
amplitude, of the channel. Alternatively the same regime can be achieved in microfluidic devices by exploiting low 
polar solvents characterized by larger Debye lengths, e.g. A ~ 1.6/im for CHB-decalin 0 Accordingly, Fig.[2]b shows 
that negatively charged tracers experience larger permeabilities for (j) = 0.5 than for (j) = 0 (compare the blue big 
(small) triangles and blue big (small) circles) while the opposite holds for positively charged tracers that experience a 
larger permeability when the bottleneck is at the middle of the channel (cj) = 0). The qualitatively different response 
of positively and negatively charged tracers to channel corrugation indicates the relevance of the geometrical coupling 
together with the electrostatic attraction or repulsion to the walls; such a dependence is encoded in the free energy 
dependence, AA. The analysis at the end of the previous section illustrates how a change in the sensitivity to spatial 
variations of the effective free energy has a strong impact in tracer diffusivity. The results obtained highlight the 
relevance of the regime where the Debye length and channel section are comparable in size. 

Since tracers move with different effective diffusivities along the channel, depending on their charge, tracer motion 
as a result of the applied density gradient will induce also a net electric current. This is an electric current induced 
by the tracers because the electrolyte is in equilibrium, and corresponds to a coupled transport effect induced only by 
geometrical variations [^ . For simplicity, we will consider that the channel is subject to the same density gradient 
for negative, positive and neural tracers, A/r+ = A^_ = A/tq. In this case, rather than computing the electric current 
itself, / = J+ — J_, it is insightful to analyze the tracer electric current relative to the reference mass current of 
neutral tracers. 


^ x+ - X- 

Jq Xo 


(24) 


which, according to Eq. (HZl, can be expressed in terms of the corresponding effective tracer permeabilities. Fig. [21c 
shows the relative electric current as a function of the position of the maximum aperture of the channel. The results 
sow that a relevant current can be induced by the tracer density gradient and that the current can change its sign 
depending on the channel geometry. The channel where the relative electrostatic current is positive correlates with 
the situations where the diffusion of positive tracers is enhanced, as could be expected. Therefore, it it is possible to 
control the magnitude and the sign of the electric current by properly tuning the position of the channel bottleneck. In 
particular, the maxima of the electric current are obtained for cj) = 0, cj) = tt, i.e. when the bottleneck is at the center 
of the channel or at channel edges and the absolute value of the current is maximized in the bottleneck case {(j) = 0). 
Hence, the shape of the channel controls the amplitude of the net electric current induced by the tracers density 
contrast imposed at channel ends. Fig. reffig:current.d shows that the amplitude of the current is modulated by the 
Debye length, A, and the maximum is achieved when is comparable to the characteristic channel section. Fig. |21d also 
shows that this current vanishes when the double layer is much thinner or wider than the channel aperture. In this 
situation the inhomogeneous electrolyte distribution in the channel is negligible. This feature highlights the fact that 
the induced tracer current emerges from the combined geometric modulation and the charge distribution across it. 
Therefore, varying the channel geometry and the electrolyte ionic strength it is possible to control both he magnitude 
and the sign of the electric current 

The modulation in channel geometry as well as the magnitude of the Debye length affects also the tracer distribution 
along the channel. Fig. [3] shows that even for neutral tracers, for which the density profile is independent of A, the 
density profile is affected by the geometrical constraint. In particular, in the case of a flat channel, the chemical 
potentials contrast imposed at the ends of the channel leads to a linear density profile inside the channel (data not 
shown). On the contrary, when the section of the channel varies, we observe a departure from the linear profile that is 
asymmetric with respect to channel geometry. In particular, for (j) = Q (cj) = 0.5) we observe an excess of accumulation 
(depletion) of particles close to the reservoir characterized by the larger chemical potential while the opposite holds 
in the vicinity of other reservoir. Interestingly, the excess of particle density is sensitive to the derivative of channel 
section; the maximum of the deviations from the linear profile is obtained when dxh(x) is maximum. 

When tracers are charged, the effect of the confinement on tracers’ density is amplified. Fig|4| shows an enhanced 
accumulation (depletion) as compared to the case of neutral tracers. In particular, positively charged tracers, (at¬ 
tracted by the negatively charged channel walls) tend to accumulate for (j) = 0, i.e. when the channel bottleneck is at 





FIG. 3: Neutral tracers density profile, normalized by the density profile in the case of a flat channel p{hi = 0), for different 
channel geometries AS' = 1.1, 2, 3.7 for squares, circles and triangles respectively characterized by, (j> = 0 (left), (p = 0.5 (right). 




x/L x/L 

FIG. 4: Charged tracers density profile, normalized by the density profile in the case of a flat channel p{hi = 0), for positive 
(red markers) and negative (blue markers) tracers and different channel geometries AS = 1.1, 2, 3.7 for squares, circles and 
triangles respectively characterized by, 0 = 0 (left), 0 = 0.5 (right). 


x = 0 while the opposite holds for negatively charged tracers. The amplitude of the excess accumulation/depletion 
of tracers has a non-trivial dependence on channel geometry. In particular, comparing the two panels of Figl4]we see 
that the density profiles for the two cases, namely 0 = 0 and 0 = 0.5 have different shapes and the amplitude of the 
modulation in the density profile is larger for 0 = 0.5. The electrolyte properties play a relevant role in determining 
tracer density profiles. As shown in Fig. |4l the maximum amplitude of the modulation in the density shows a non- 
monotonous behavior with A and its maximum is obtained for khf) ~ 10, i.e. when the amplitude of the bottlenecks, 
k{hQ — hi) ~ 0.5 is comparable to the Debye length. 


V Conclusions 

We have studied the impact of the channel shape on the permeability of charged and neutral tracers. We have 
identified a novel mechanism that controls tracer permeability due to variations in the channel geometry. In particular, 
we have studied the geometrical dependence of the channel permeability for neutral and charged tracers, when the 
system is driven by means of a chemical potential contrast imposed by controlling tracers concentration at channel 
ends. In order to keep analytical insight, we have assumed that the channel amplitude varies smoothly, dxh{x), 
along its longitudinal axis, x. In such a regime, a lubrication approximation is reliable and allows for a significative 
simplification of the analysis. Moreover, in such a regime it is possible to exploit the Fick-Jacobs approximation 
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that allows us to factorize the distribution probability of tracers, hence reducing the problem to a \D problem whose 
solution can be accessed analytically. 

We have exploited such a framework to characterize the permeability of a corrugated channel for the transport 
of positive and negative charged and neutral tracers. For these systems, we have found that the relative magnitude 
of the permeabilities of both positively and negatively charged tracers as well as neutral tracers can be tuned by 
the geometry of the channel, leading to a variety of scenarios where positively (negatively) charged tracers, can 
experience larrar or smaller permeabilities than neutral ones. As it has already been shown in the case of electric 
driving forces [^, the geometrically-induced control on tracers current is sensitive to electrolyte properties captured 
by the Debye length, A. The sensitivity of tracer diffusion to channel corrugation is maximum when A is comparable to 
the characteristic channel section. Such a regime kh^ ~ 1 can be obtained in nanofluidic devices for which the Debye 
length, that is typically between 1 and 10 nm, matches the typical amplitude, /iq of the channel, or in microfluidic 
devices by exploiting low polar solvents characterized by larger Debye lengths such as A ~ l.h/iw [l^. In contrast, 
for very wide and narrow double layers, i.e. for kho ^ 1 or kho <C 1, the impact of channel corrugation vanishes 
asymptotically. The interplay between channel corrugation and electrolyte structuring across it leads to a new cross 
transport mechanism. Specifically, we have found that when the Debye length is comparable to the channel section, a 
net electric current develops when a chemical potential gradient of positive and negatively charged tracers acts at the 
channel ends. This corrugation-induced cross transport effect can be of significance in the electric transport through 
ionic channels and deserves to be analyzed in detail. 
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Appendix A 


In order to keep the argument of the logarithm in Eq. [3 dimensionless, we have introduced a prefactor, l/2hoLj;. 
In the following we show that the choice of this prefactor is arbitrary since it will not affect either the probability 
distribution pa or its flux. The starting point are Eos. I2l6l7l that we report here: 

dtPa{x,y,z,t) = D/3V ■ {Pa{x,y,z,t)VUa{x,y,z)) + DV'^Pa{x,y,z,t) 


with the ansatz: 


^-I3qceip{x,y,z) 

Pa{x,y,Z,t) = Pa{x,t) - -l 3 A^{x) - 


PAaix) = — In 


1 


f>Lzl2 ph{x) 


“^hoLz J-L^/2 J-h(x) 




We can rewrite the last equation as j3Aa{x) = —In 


2hQLz 


- In 




(A.25) 

(A.26) 

(A.27) 

From 


the last expression we can calculate e = e P^°e where Aq = In 


^HqL z 


and Ai(x) = 


In 


S-LJ 2 f-it) e . Substituting the last expressions in Eq.l and integrating in y we get: 


-Lz /2 J —h{x 

r Q-Pqceit{x,y) 

P / dy pAq,-PA,ix) = Dpdx 


Paix,t) 


g-Pqceip{x,y) 

g-PAog-PAi(x) 


d^U + Ddlpa{x,t) 


^-Pqc,etj}{x,y) 

g—PAoQ-PAi{x) 


dy 


(A.28) 


where in the last step we have neglected the terms like dy J[...]dy since they provide null contributions. From the last 
equation it is clear that we can get rid of the term e~P^° since, being independent on both, x and y, it can be taken 
out of the integrals and differentiation operators. Recalling that J dy *^ = 1 and performing the integral in 

the last expression we get: 


Pa{x,t) = dxD 


^ , .dAiix) „ , , 

dPa {X, t) --- - + dxPa (X, t) 


dx 


Since dxAi{x) = dxAa{x) we can rewrite the last expression as 

Pa{x,t) = dxD 


I3paix,t)^^^^ + dxPa{x,t) 


dx 


(A.29) 


(A.30) 
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showing that the normalization does not affect the evolution of the tracer probability along the channel. 
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